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Abstract 

The nuclear (re; t)— reaction of the three-electron 6 Li atom with thermal/slow neutrons is con- 
sidered. An effective method has been developed for determining the probabilities of formation of 
various atoms and ions in different bound states. We discuss a number of fundamental questions 
directly related to numerical computations of the final state atomic probabilities. A few appro- 
priate variational expansions for atomic wave functions of the incident lithium atom and final 
helium atom and/or tritium negatively charged ion are discussed. It appears that the final 4 He 
atom arising during the nuclear (n, 6 Li; 4 He, t)-reaction in the three-electron Li atom can also be 
created in its triplet states. The formation of the quasi-stable three-electron e% during the nuclear 
(n; t)— reaction at the Li atom is briefly discussed. Bremsstrahlung emitted by atomic electrons 
accelerated by the rapidly moving nuclear fragments from this reaction is analyzed. The frequency 
spectrum of the emitted radiation is investigated. 
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I. 



INTRODUCTION 




6 Li + n = 4 He + 3 H + 4.785 MeV 



(1) 



where the notations 4 He and 3 H stand for the helium nucleus (also the called the a— particle, 
or a, for short) and tritium (or t) nucleus. The cross-section o of this nuclear reaction for 
thermal neutrons with E n « is very large a max ~ 960 • 10~ 24 cm 2 (or 960 barn) 
The velocities of the two nuclear fragments formed in the reaction, Eq.flT]), with thermal 
neutrons are Vt ~ 6.03986 a.u. and v a ~ 4.52989 a.u. for the tritium nucleus and a— particle, 
respectively. In this study all particle velocities are given in atomic units, where H — 1, m e = 
1, e = 1 and the unit of atomic velocity is v e = ac ~ 2.1882661 • 10 8 cm • sec -1 . Here 
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and everywhere below c is the speed of light and a = j- is the dimensionless fine structure 
constant. This 'atomic velocity' v e is the velocity of the Is— electron in the hydrogen atom 
with the infinitely heavy nucleus °°H. It is clear that in atomic units v e = 1. 

In reality, at normal conditions the nuclear reaction, Eq.(l), occurs in the three-electron 
6 Li atom. Therefore, the electron density of the original Li atom is re-distributed somehow 
between the two new nuclei created in the reaction, Eq.([T]). The newly created atomic nuclei 
move rapidly with velocities which significantly exceed average velocities of atomic electrons 
in the incident Li atom. The re-distribution of electron density of the incident atom between 
two rapidly moving nuclei leads to the three following processes: (1) formation of the final 
atoms and/or ions in a variety of different excited states, (2) emission of breaking radiation 
(or bremsstrahlung) by the accelerated electrons, and (3) formation of the spatially isolated 
quasi-stable three-electron system e% . The first goal of this study is to evaluate the final state 
(atomic) probabilities for the reaction Eq.([T]). Briefly, we want to determine approximate 
probabilities to detect the helium and tritium atoms, and He + and 3 H~ ions in their bound 
states. Note that all newly created atomic fragments from the reaction, Eq.((T]), move rapidly 
even in the case of thermal and slow neutrons. The accurate theoretical prediction of the final 
state probabilities for rapidly moving atomic fragments is not trivial. Moreover, addressing 
this problem reveals a large number of unanswered fundamental questions. In this study we 
have discovered, among other things, that the final 4 He atom can be formed not only in the 
singlet spin states, but in a number of triplet states also. Another interesting possibility is the 
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formation of the quasi-stable three-electron system during the nuclear reaction, Eq.flT]). 
Such system e% contains no positive (heavy) particles and has a number of unique properties. 
In this study we also want to analyze bremsstrahlung emitted after the nuclear reaction, 
Eq.flTJ), in the three-electron Li atom. It is shown below that the frequency spectrum of such 
radiation essentially coincides with the well known bremsstrahlung spectrum. 

II. SUDDEN APPROXIMATION 

As we have mentioned above the rate of the nuclear reaction, Eq.flTJ, is much faster than 
the rates of atomic transitions in the incident Li atom. Moreover, the two new nuclei ( 4 He 
and 3 H) arising in this reaction move with the velocities which are significantly larger than 
the corresponding electron velocities in the incident Li atom. This means that to determine 
the final state probabilities (i.e. probabilities to find newly formed atomic systems in some 
final states) we can apply the so-called sudden approximation 0], j^. Note also that 
after the nuclear reaction, Eq.flT]), in the three-electron Li atom we have to deal with three 
new sub-systems. Two of these systems are the positively charged 4 He and 3 H nuclei rapidly 
moving from the reaction area (fast sub-systems). The third (slow) sub-system includes 
three electrons remaining in the reaction area, i.e. the three-electron (or tri-electron) . 
These three electrons are left from the incident Li atom and their wave functions still keep 
information about the original atom which disappeared during the nuclear reaction, Eq.flT]). 
Interactions between these three sub-systems produce all phenomena mentioned above, i.e. 
formation of the different atomic species in a variety of bound (and unbound) states and 
emission of breaking radiation. 

Suppose we want to determine the probability that the rapidly moving 4 He nucleus takes 
one of the three atomic (Li) electrons remaining unbound after the nuclear reaction, Eq.flT]). 
In other words, the final atomic system is the positively charged helium ion 4 He + . For an 
atomic system containing fast and slow sub-systems the corresponding probability is written 
in the form (see, e.g., Q]) 

Pif =| J Pu(r u )exp(tV a ■ ri 4 )p H c+(ri4)d 3 ri4 | (2) 

where V Q is the velocity of the He + ion after the reaction, Eq.flT]) and Pu{ r u) an d PHe+( r i4) 
are the one-electron density matrixes of the Li atom and He + ion, repectively. These density 
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matrices are simply related with the corresponding wave functions of the bound states, e.g., 
for the one-electron density matrix of the Li atom 

Phi( T u) = ^J J J ^Li( X 14,X24,X 3 4)^ L i(xi4,X24,X3 4 )(iSi(iX24(iX34 (3) 

where = (sj,rj) is the complete set of spin-spatial coordinates of the i-th electron. In 
Eqs.((2]) - ()3]) and everywhere below the nucleus in the three-electron Li atom is designated 
by the index 4, while three electrons are denoted by the indexes 1, 2 and 3. This definition 
of the one-electron density matrix, Eq.(j3])), and notations used in this equation correspond 
to the definition given in Eqs.(4.1.5) - (4.1.6) of jfj]. The one-electron density matrix for the 
one-elecron He + ion is 

PHe+0l4) = P ( nl+( r u) = J $Hc+( X 14)$Hc+(xi 4 )d.Si (4) 

It contains only one integral over spin variables of the electron 1. 

If the tritium atom is formed during the nuclear reaction, Eq.flT]) in the three-electron Li 
atom, then we can write the following formula for the final state probability Pif 

Pif = J Pu ( r u) exp(W t ■ ri 4 )p T (ri 4 )<i 3 ri 4 (5) 

where Pt(i"i4) is the one-electron density matrix of the tritium atom, and V$ is the velocity 
of the tritium nucleus after the nuclear reaction, Eq.([T]). Here the notations T and t stand 
for the tritium and and tritium nucleus, respectively, while V t is the velocity of the tritium 
nucleus after the nuclear reaction, Eq.flT]). 

For the helium atom product (two bound electrons) the final state probability Pif is 
written in the form 



if 



pi- (r M , r 24 ) exp(tV a ■ r u + iV a ■ r 24 )pHe (r 14 , ^c^r^cf^ (6) 



where V Q is the final velocity of the a— particle after the nuclear reaction, Eq.([T]). The 
two-electron density matrix of the three-electron Li atom is defined as follows 

PL i ) (ri4,r 24 ) = Qj J J ^Li( X 14,X 2 4,X3 4 )^ Li (xi 4 ,X24,X3 4 )c?SiC?S2rfX3 4 (7) 

where we have used the factor iV(iV — 1) in the front of the density matrices In some 
cases, however, it is better to use a different definition of the density matrices with the unity 
factor. The two-electron density matrix of the He atom (bound states) takes the form 

pHe(ri4, r 24 ) = PhI ( r i4, r 24 ) = / $H e ( x i4, r 2 4)$He(xi4, T 2 A)ds 1 ds 2 (8) 



In the case of nuclear reaction, Eq.([T]), in the three electron Li atom we need only one- 
and two-electron density matrices, since the negatively charged He~ is unstable. Note that 
all these integrals mentioned in the expressions for the probabilities, Eq.(j2]), Eqs.(j5]) - flSJ) 
and density matrices, Eqs.OH]) - (jlj) and (jSJ), include the wave functions of the bound few- 
electron atoms and ions. These wave functions must be determined from accurate atomic 
computations, whose construction is discussed in the next Section. 



III. VARIATIONAL WAVE FUNCTION OF THE LITHIUM ATOM 



To determine the final state probabilities we need to construct the accurate variational 
wave functions for the incident and final atomic systems involved in the reaction, Eq.flT]). 
In our earlier paper |8j] we have calculated the final state probabilities to form one-electron 
atoms and ions in those cases when exothermic nuclear (n;t)— and (n; a)— reactions occur 
in one-electron atoms/ions, i.e., both incident and final atomic systems contain only one 
bound electron. In this study we deal with the actual three-electron wave function of the 
incident 6 Li atom and two-electron wave function of the final 4 He atom. Also, we consider the 
formation of the one-electron He + ion, tritium atom 3 H (or T) and two-electron negatively 
charged tritium ion 3 H~ (or T~). 

First, let us discuss the construction of three-electron variational wave function for the 
Li atom. Without loss of generality, below we restrict ourselves to the consideration of the 
ground 2 S(L = 0)— state of the Li atom. As is well known (see, e.g., js], 10]) the accurate 
variational wave function of the ground (doublet) 2 S(L = 0)— state of the Li atom is written 
in the following general form 

\l/(Li) L= o = ip L=0 (A] {r»j})(a/3a - patat) + 4> L=0 (B] {r ij }){2aaf3 - (5aa - afta) (9) 

where if>L=o(A; {r^}) and (pL=o(B; {r^}) are the two independent radial parts (= spatial 
parts) of the total wave function. Everywhere below in this study, we shall assume that all 
mentioned wave functions have unit norms. The notations a and /3 in Eq. stand for the 



one-electron spin-up and spin-down functions, respectively (see, e.g., III]). The notations 
A and B in Eq.Q mean that the two sets of non- linear parameters associated with the 
radial functions if) and (f> can be optimized independently. Note that each of the radial basis 
functions in Eq.Q explicitly depends upon all six interparticle (or relative) coordinates 
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r i2, ^13, r 23, r i4, "^24, "^34, where the indexes 1, 2, 3 stand for the three electrons, while index 4 
means the nucleus. 

The actual atomic wave function in an atomic system must be completely antisymmetric 
with respect to all electron spin-spatial variables. For a three-electron wave function this 
requirement is written in the form ^4.123^(1,2,3) = —^(1,2,3), where the wave function 
\I/ is given by Eq.(j9]) and A e is the three-particle (or three-electron) antisymmetrizer A e = 
e—Pi2—Pi3—p23+Pi23+Pi32- Here e is the identity permutation, while Py is the permutation 
of the i-th and j-th particles. Analogously, the operator P^ is the permutation of the i-th, 
j-th and k-th particles. 

Suppose that the three-electron wave function of the Li atom has been constructed in 
the form of Eq.Q. By applying the antisymmetrizer *4i23 to the first part of the total wave 
function, Eq.([9]), one finds 

-4-123 ipL=o{A; {rij})(aj3a — j3aa)] = (etp)(a(3a — (3aa) + (Pi2ip)(ocfioc — f3aa) 
— (Pi 3 ip)(a(3a — aa(3) — (P 23 ?/>) (aa/3 — (3aa) + (Pi23ip)(aa(3 — a(5a) 

+ (Pi324>)(f3aa - aa(3) (10) 

where the notations (Pijip) and (Pijkip) mean the permutation operators which act on the 
coordinate wave functions only. Analogously, for the second part of the total wave function 



one 



finds 



-4-123 4>l=o(B; {rij})(2aa(3 — (3aa — a(3a)] = (£(j))(2aa(3 — (3aa — aj3a) 
— (Pi20)(2aa/3 — (5aa — afia) — (Pi 3 0)(2/9aa — aa(5 — a(3a) 
— (P230)(2a/3a — aa(3 — f3aa) + (Pi230)(2/3aa — a(3a — aa(5) 

+ (Pi320)(2a/3a — aa(3 — f3aa) (11) 

where the notations {Pij<p) and (Pijk<fr) mean the permutations of the spatial coordinates in 
the (pL =0 (B; {r^}) radial function, Eq.([H]). 

Now, by using the formulas, Eqs.f llOp and ffTT]) . we can obtain the formulas which can 
be used in computations of the final state probabilities in the case of nuclear reaction, 
Eq.([l]), in the three-electron Li atom. For instance, if the final wave function has the same 
spin-symmetry, i.e. it is written in the form 

^fi — ipfi(ri, r 2 , r 3 )(a/3a - f3aa) + 0/i(r x , r 2 , r 3 )(2aa/3 - (3aa - a {3 a) (12) 



then the final state probabilities are determned with the use of the following formulas 

P^ = (^(n, r 2 , r 3 ) | ^(2e + 2P 12 - P 13 - As - P 123 - P 132 )M^ { nj })) (13) 

P Hl =(<t>fi(T 1 ,T 2 ,T 3 ) | ~(Pl3- A3 + A23- A 32 )^Li(A;K})) (14) 

P^ = ftMri,r 2 ,r 3 ) | J (As - A3 + A23 - A32WOB; {r y })> (15) 



P 00 = (0 ft ( ri , r 2 , r 3 ) I ^(2e - 2A2 + A3 + A3 - A 23 ^ Pisa) MB] {%})) (16) 



12| which correspond to the 



Note that these formulas coinside with the known formulas 
doublet — > doublet transition in the three-electron atomic systems. Briefly, this means 
that the both incident and final atomic states contain three-electrons in the doublet spin 
configuration (the total electron spin equals |). 

By considering all possible spin configurations for the final atomic state one finds that 
other spin configurations are possible. For instance, if reaction Eq.flTJ leads to the formation 
of the two neutral atoms ( 4 He and 3 H atoms), then the helium atom can be detected either 
in the singlet state, or in the triplet state. The spin function of the singlet state is Xs — 
4g(a;/3 — /?«), while the wave function of the triplet state can be chosen (in our case) in 
the form Xt — olol. It is clear that these two wave functions have unit norms and they are 
orthogonal to each other. The probability of formation of the 4 He atom in its singlet spin 
state is determined by the fourmulas, Eqs. ffTB"]) - ffT6]) . The tritium atom in this case will be 
formed with the /3— electron spin function. 

If the final state of the 4 He atom is the triplet state (tr), then the formulas for the final 
state probabilities take the form 

Ptr^= (<Mri,r 2 ,r 3 ) | ^(A 3 - As + A23 - A 32 )^Li(^ {r y })) (17) 

(M^ r 2,r 3 ) I ^(2e-2P l2 + P 13 + P 23 -P 123 -P 132 )MB;{r ij })) (18) 

These formulas indicate clearly that the probability to find the 4 He atom arising during 
the nuclear reaction, Eq.flT]), in its triplet spin state(s) is not zero. In all earllier studies the 
transitions to the final atomic states with different spin states was ignored (never considered). 



IV. FORMATION OF THE TRI-ELECTRON 



As mentioned above the two nuclear fragments formed in the reaction, Eq.flTJ), move 
with very large velocities. They leave the reaction zone very quickly for a time which 
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is significantly shorter than a typical atomic time r a = ~ 2.42 ■ 10~ 17 sec. The three 
negatively charged electrons e~ remain in the reaction zone substantially longer. The arising 
quasi-stable system of three electrons is called a three-electron, or tri-electron. It is clear 
that the tri-electron is not a stable system and after some time it will be destroyed by 
the Coulomb repulsion between electrons. However, in some cases such a destruction takes 
a relatively long time which is sufficient to detect it experimentally and measure some of 
its properties. A combination of the stabilizing electric and magnetic fields can be used to 
increase the life-time of a tri-electron. Some of the properties of the e% system can be studied 
by using modern methods of the femto- and attosecond physics. Experimental interest to 
study such new systems cannot be overestimated. 

The 'quasi- atomic' tri-electron system has a large number of unique properties and can 
be considered as a light quasi-atom which has no central (i.e. heavy), positively charged 
nucleus. In the case of reaction, Eq.flTJ, the spin-spatial symmetry of the three-electron 
wave function exactly coincides with the corresponding symmetry of the wave function of 
the ground 2 S'-state of the Li atom. 



V. CALCULATIONS 



In this study we have developed an approximate procedure to construct the bound state 
wave functions of the three-electron atoms and ions. This procedure can be used to perform 
numerical evaluations of the final state probabilities in the case of the nuclear reaction, 
Eq.([ID, and for the /^"-decaying isotopes of the three-electron atoms [3]. In this approach 
the trial wave function is constructed as the sum of various terms and each of these terms 
contains the products of the electron-nucleus functions. None of the three electron-electron 
coordinates r 32 , r 31 , r 2 i is included in such trial wave functions. This simplifies drastically 
the following Fourier transforms of such wave functions. For the ground state (the doublet 
2 S(L = 0)— state) of the Li atom the radial wave function ipL =0 (A; {r^}) is chosen in the 
following form: 

ipL=o(ru, r 2 4, r 34 , 0, 0, 0) = c 'fc r u r 24 »m 3{k 'exp(-a k r u - (3^24 ~ I^m) (19) 

fc=i 

N 3 

= Yl C' fc r™ l(fc) r™ 2(A:) r™ 3{fc) exp(-a fc ri - (3 k r 2 - 7 fe r 3 ) 
k=i 
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where C/~ are the linear (or variational) coefficients, while mi(k), rri2(k) and m^{k) are the 
three integer (non-negative) parameters, which are, in fact, the powers of the three electron- 
nucleus coordinates r i4 = n (i = 1, 2, 3). The real (and non-negative) parameters ak,Pk,1k 
are the 3N S varied parameters of the variational expansion, Eq. (IT9"|) . Below, we shall assume 
that the trial wave function Eq. fl!9j) has a unit norm. Furthermore, in all calculations 
performed for this study only one spin function = a P a ~ 0oca) is used. The wave 

function, Eq. ffT9]) . must be properly anti-symmetrized upon all spin-spatial coordinates of 
the three electrons. 

The principal question for the wave function, Eq. ffTDl . is related to its overall accuracy. 
If (and only if) such accuracy is relatively high, then such a wave function, Eq. (1191) . can be 
used in actual computations of the probability amplitudes. In this study we have constructed 
the 23-term variational wave function shown in Table I. This wave function is represented 
in the form of Eq. (jT9l) and contains no electron-electron coordinates. All sixty nine (69 = 
3 x 23) non-linear parameters 7fc (k = 1, 2 . . . 23) in this wave function have been 

optimized carefully in a series of bound state computations performed for the ground state 
of the Li atom. Finally, the total energy E of the ground 2 S— state of the °°Li atom obtained 
with this 23-term trial wave function, Eq. (fT9l) . was -7.44859276608 a.u. This energy is close 
to the exact total energy of the ground state of the °°Li atom. It indicates a good overall 
quality of our approximate wave function with 23 terms which does not include any of the 
electron-electron coordinates ?"i2, ri3, r23. This wave function is used in the computations of 
the final state probabilities (see below) for the nuclear reaction, Eq.Q, in three-electron Li 
atom. 

Note also that in atomic physics based on the Hartree-Fock and hydrogenic approxima- 
tions the ground state in the Li atom is designated as the 2 2 S— state, while in the clas- 
sification scheme developed in highly accurate computations the same state is designated 
as the 1 2 S— state because it corresponds to the first eigenvalue of the Hamiltonian matrix 
H. This classification scheme is very convenient for truelly correlated few-electron wave 
functions which represent situation where there are no good hydrogenic quantum numbers. 
Nevertheless, to avoid confllicts between these two classification schemes in this study we 
follow the system of notation used earlier by Larsson [9] which designated this state of the 
Li atom as the 'ground 2 S , -state'. 

The advantage of the basis functions, Eq. (1191) . for numerical calculations of the final state 
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probabilities for the reaction Eq.flT]) is two-fold: (1) these trial wave functions allow one to 
obtain relatively accurate approximations for atomic three-electron wave functions; and (2) 
all required Fourier transforms are easily perfomed for such trial wave functions since they 
do not include electron-electron coordinates. Note also, that approximate two-electron wave 
functions of the He atom and/or ion can be constructed in analogous form which contain 
only electron-nucleus coordinates, e.g., 

i N s 

^=o(ri, r 2 , 0) = -=(1 + P u ) £ C k r? lW rF {k) e*P{-<* k ri - to) (20) 
V 2 k=1 

Some probabilities (in %) of observing selected electron final states in the helium ion 
4 He + and tritium atom 3 H arising in the exothermic nuclear reaction Eq.flT]) of the three- 
electron lithium-6 atom can be found in Table II. In calculations included in Table II we 
have assumed that the incident Li atom was in its ground 2 S— state. Also, we assume that 
the nuclear (n; t)— reaction, Eq.flT]), proceeds is produced by thermal and/or slow neutrons. 



VI. BREMSSTRAHLUNG 



As we have mentioned above the two positively charged atomic nuclei created in the 
nuclear reaction, Eq.flT]), with thermal neutrons move with substantial velocities in the two 
opposite directions. In general, these two positively charged nuclei will produce acceleration 
of electrons which still remain in the reaction zone. In turn, the accelerated electrons will 
generate emission of radiation. It appears that the spectrum of this radiation essentially 
coincides with the known spectrum of bremsstrahlung (or breaking radiation). In this Section 
we discuss the basic properties of radiation emitted during the nuclear reaction, Eq.flT]), in 
three-electron 6 Li atom. 

Note that the reaction, Eq.flT]), is a fission-type reaction in a three-electron atomic system. 
The original nuclues decays into two electrically charged fragments which are rapidly moving 
away from each other. The general theorv of radiation emitted during such processes was 
developed in our earlier studies (see, e.g., [y| and references therein). Atomic electrons 
from outer electronic shells become free during this fission-type nuclear reaction in an atomic 
nucleus. These free electrons interact with the rapidly moving nuclear fragments and such 
an interaction can produce breaking radiation or bremsstrahlung. Here we want to derive 
the explicit formulas for this radiation and investigate its spectrum. The electric charges 
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and velocities of these fission fragments are Qie, Q 2 e and Vi, V2, respectively. To simplify 
all formulas below, we shall assume that the both fission fragments move along the Z— axis. 
Furthermore, in this study we restrict ourselves to the consideration of the non-relativistic 
processes only (i.e. the two velocities Vi, V% are significantly less than the speed of light c in 
vacuum). The case of arbitrary velocities is discussed in |14j |. 

The electron acceleration is directly related to the second time-derivative of the dipole 
moment, i.e. d = er, where d is is the dipole moment. In the case of fission-type reaction 
in atomic systems the explicit formula for electron's acceleration is 

^ " 2 -- Qie*/ 1 x , Q 2 e 2 



m, 



(21) 



where R\ = J x 2 + y 2 + (z — V\t) 2 and R 2 



(z + V 2 t) 2 . The notations V\ and 



Vi stand for the velocities of the two fission fragments. The second time-derivative of the 
dipole moment is 



Qie 



v(-L) + 9*l v{ ±. 



Qie 3 Ri Q 2 e 3 R 2 



R 3 , 



(22) 



m e ^tix J m e x ti 2 ' m e R\ m e 

where Ri = (x, y, z — V\t) and R2 = (x, y, z + V^t) are the three-dimensional vectors. The 
intensity of the non-relativistic bremsstrahlung from a fission-type reaction in a one-electron 
atomic system is 



dl 



(d x nfdQ = (—) 



e 2 N2QV r QiRi Q 2 R 2 



sin 9dQ 



(23) 



4ttc 3V ^"" / ^mj 4vrc 3 Lg Rf Q R\- 

where 9 is the angle between the vector d and vector n which designates the direction of 
propagation of radiation. The notation Q in the last equation stands for an arbitrary electric 
charge. This value can be considered as an output parameter. In particular, one can choose 
Q — Qii or Q = Qi. For a iV e — electron atomic system the last formula takes the form 

„2 . „ i\T r>2„2 



dl 



(-) 



2N e Q 2 e 2 rQi Ri Q2R2 
Q R\ + Q Rl 



Aire 3 



sin OdVt 



This formula can be re-written in the form 



dl ( c \ 
dQ ^ m e ' 



e 2 ^N e Q 2 e 2 r Q| J_ Q|J_ Q X Q 2 ^ ■ R 2 
Q 2 Rf Q 2 Ri 



Aire 3 



Q 2 V #3^,3 

for the differential cross-section. Our derivation of these formulas is based on the fact that 
radiation emitted by different post-atomic electrons is non-coherent. Indeed, the three free 
electrons arising during the reaction, Eq.(TfJ), are independent quantum systems and they 



sin 2 9 



(24) 



(25) 
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move as 'random' particles. It can be shown that the 'randomization' of the phases of 
the electron wave functions is directly related to the electron-electron repulsion which was 
ignored in Eqs. ()2ip - (I2"5j) . At certain conditions bremsstrahlung from fission- type process 
in atomic systems can be coherent. For instance, if two fission fragments move with very 
large velocities, while all post-atomic electrons almost do not move at the beginning of the 
process. In this case, one needs to introduce an additional factor N e in the last formula. 

Let us discuss the spectrum of the emitted radiation. As follows from the formula, 
Eq.( )25|) . the intensity of bremsstrahlung from fission-type processes in atomic systems 
rapidly decreases with the time I ~ t~ A after the nuclear reaction (t — 0). Such a behaviour 
is typical for all finite-time process/reactions with accelerated electrons. The spectral res- 
olution R(u) (or spectral function) of the intensity of dipole radiation is written in the 
form 

. 4 , - , 9 <ko 4a; 4 . , , 9 duo , . 

dR{u) = — dj 2 — = — <L 2 — 26 
3cr 2tt Scr lit 

where d w is the Fourier component of the dipole moment d defined above. As follows from 
Eq. fl2"61) one needs to find the Fourier components of the vector of the dipole moment d^, 
Eq. (l21j) . Finally, the problem is reduced to the calculation of the two following Fourier 
transformations 

h(uj- -,V, cos 77) = / — ±2(-cosr/) -t + t 2 2 exp(tut)dt (27) 
V Jo l V z V J 



and 



q z'+oo r q2 q 3 

I 2 (uj;—,V,cosr})= —r±2(-cosr])-t + t 2 2 texp(tut)dt (28) 

V Jo l V V J 



where a is the 'effective' radius of the original electron shell and r\ is the angle between 
the electron acceleration and Z— axis (or the line along which the two nuclei are moving). 
Note that the spectrum of the emitted radiation depends upon cos 77, as for all fission-type 
reactions and processes. The lower limit in these formulas is zero, but, in reality, we cannot 
use times t which are shorter than r = -. Formally, this means that the lower limits in 
Eqs. (I27j) and ff28l can slightly be changed and this can be used to simplify the explicit 
formulas for the corresponding Fourier components. The actual upper limits in Eqs. (l27|) 
and (1281) are also finite. For instance, it is possible to obtain a very good approximation for 
the integrals, Eqs.( l2T|) and (12"B"1) . by using the upper limit In actual situations one can 
use the formula exp (ix) = cos x + 1 sin x and then calculate all arising integrals numerically, 
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or by using analytical formulas from the Tables of Fourier sine/ cosine transformations (see, 
e.g., (J). 

To produce the explicit formula for the spectral functions h(co; V, cos 77) and 
I2(oj; y, V, cosrj) we note the following relation between these two functions: 

& d (x 

h(u; y, V, cosr/) = -i—h(u; — , V, cost?) (29) 

This formula substantially simplifies the computation of the Fourier transformations. More- 
over, a very good approximation for the y, V, cosn) function is: 

h(u; V.cosij) = ex P (-za;-,cos^ [-KiM + -ln( T -|^-) - I (1 _^ )2 j (30) 

where b = aS y V and Ki(x) is the MacDonald function (see, e.g., jl6[ ) . The formula, Eq.( 13"U|) . 
describes the spectrum of radiation emitted by one electron accelerated by rapidly moving 
atomic fragments from the nuclear (n, 6 Li; 4 He, £)— reaction. The actual spectra emitted 
during the nuclear (n, 6 Li; 4 He, t)— reaction in the three-electron Li atom can be different 
from the spectrum predicted by the formula Eq. fl30|) since we have ignored all contributions 
from the electron-electron repulsion(s). It is interesting to study such radiation in detail. 
In general, the intensity of bremsstrahlung from the reaction, Eq.([T]), can be amplified by 
accelarating the incident 6 Li atom or 6 Li + ion to very large velocities. 



VII. CONCLUSION 



We have considered the nuclear reaction, Eq.([T]), with the thermal/slow neutrons in the 
three-electron Li atom. It is shown that a number of different atoms and ions can be formed 
during this reaction. This also includes the formation of the electrically charged tri-electron 
e% ion which is a quasi-stable negatively charged ion. It has the same (doublet) symmetry 
as the wave function of the incident Li atom. We have developed an effective method for 
accurate evaluation of the final state probabilities, i.e. probabilities to form different atomic 
species in a variety of bound states. This method is based on the use of special three- and 
two-electron wave functions of the Li and He atoms, respectively, whose explicit form allows 
one to obtain accurate numerical results for all bound state properties of the incident and 
final atoms/ions. On the other hand, for the same wave functions it is easy to perform 
all complete and/or partial Fourier transformations needed in calculations of the final state 
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probabilities. The spectrum of the emitted radiation has been investigated. It is shown that 
such a spectrum essentially coincides with the known spectrum of bremsstrahlung emitted 
by the accelerated electrons. 

Our procedure can now be used for the more complicated nuclear (n; a)— reaction in the 
five-electron B-atom (see Appendix). The (n, 10 B; 7 Li, a)— reaction on this atom is of great 
interest in a number of applications, e.g., in medical physics (BNCT). The Li-atom/ion 
and He-atom/ion which are formed during this reaction may contain up to three and two 
electrons, respectively. These atomic fragments move rapidly, with the velocities vu m 
2.40896 cm ■ sec~ l and v a ~ 4.21568 cm ■ sec -1 in the case of slow neutrons. The sudden 
approximation can certainly be applied to the He-atom and He-like ions. However, this 
approximation cannot be used for internal electrons (or l 2 s— electrons) of the Li atom/ion, 
since the velocities of these two electrons are comparable with the final velocity of the 7 Li 
nucleus. It is very likely that the probability to observe the 7 Li 2+ ions after Eq. (l3Tj) will be 
relatively small < 5 - 10 %, while analogous probability for the 7 Li + ions will be relatively 
large (> 50 %). 

Appendix 

In this study our analysis was restricted to the consideration of the nuclear 
(n; t)— reaction, Eq.flTJ, in the three-electron Li atom. Moreover, such a reaction is consid- 
ered for thermal and slow neutrons only. In reality, the nuclear reaction of the 6 Li nucleus, 
Eq.flTJ, proceeds with neutrons of all energies and the energy released increases almost lin- 
early with the energy of the incident neutron. For fast neutrons with E n > 1 MeV the 
energy released is substantially different from the value quoted in Eq.flTJ and the velocities 
of the atomic fragments increase correspondingly. Investigation of the final state (atomic) 
probabilities and bremsstrahlung emitted during such processes with fast neutrons is of in- 
terest in some applications. Note that the cross-section of this reaction has a large resonance 
(maximum) a ~ 4.5 barn at E n « 240 — 270 keV, but it is relatively large for neutrons of all 
energies E n < 0.8 MeV. This makes the reaction, Eq.flTJ, extremely important in the ther- 
monuclear ignition and following propagation of the thermonuclear burning wave in highly 
compressed (p > 100 g ■ cm~ 3 ) 6 LiD deuteride which is routinely used as a fuel in modern 
Lermonuclear explosive devices (see, e.g., [l?]], jl8]). The analogous (n, 3 He;p, t)-reaction 
also plays a very important role in such processes. 
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In general, the nuclear (n, t)-reactions of the 6 Li and 3 He nuclei with neutrons allow 
one to reduce drastically the overall bremsstrahlung loss from the hot combustion zone and 
increase the tritium/deuterium ratio which is crucially important to start new thermonuclear 
(d, t) —reactions. Briefly, by using the 6 LiD deuteride in modern thermonuclear explosive 
devices we can reduce the required compressions to relatively small values. In many cases 



such compressions are dozens times smaller (usually, in 25 - 40 times smaller 18]) than 
compressions required for any other (solid) thermonuclear fuel, e.g., 7 LiD deuteride. On 
the other hand, by using compressions which are provided by a standard 'primary' nuclear 
charge one can create extremely compact thermonuclear explosive devices based on 6 LiD 
deuteride. The idea to use pure 6 LiD deuteride in thermonuclear explosive devices was 
originally proposed by V.L. Ginzburg in 1948-1949 (see discussion and references in 

Some other (n;t)—, (n;p)— and (n; a)— reactions are often used in different applications. 
For instance, the nuclear reaction of 10 B nuclei with slow neutrons 

10 B + n = 7 Li + 4 He + 2.791 MeV (31) 

is extensively used in the boron neutron capture therapy (BNCT, for short), or boron 
neutron-capture synovectomy [3] - {23], to treat different forms of cancer, including brain 
cancer. The fast a— particle produced in the reaction, Eq.( l3ip . kills (or at least 'badly 
damages') one cancer cell before it finally stops. The modern applications of this reaction 
to cancer treatment are based on the use of molecules which contain a large number of 
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B-atoms, e.g., the IN 



molecular clusters 



2l| 



a 3 



Bon H17NH9 CH9CH9 NHol molecule, other similar molecules, and 



20 11171M12 yiio ^iio, in±±2J 
22] (see also 23( and references therein). In this case the overall 
energy release from the reaction, Eq.( j3TT) . in one cancer cell can be extremely large. Corre- 
spondingly, the local temperature in the whole cell suddenly increases to very large values 
and this kills the incident cancer cell with almost 100 % probability. Note that the tritium 
nucleus does not form in the nuclear reaction Eq. pTj) which means that it is safe to ini- 
tiate this reaction inside of a human body. By studying the nuclear reaction, Eq.([T]), in 
few-electron atoms and ions we want to develop a number of reliable theoretical methods 
and numerical procedures which can be later used in applications to the analogous reaction, 
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TABLE I: An example of the trial, three-electron wave function constructed with the use of N = 23 
semi-exponential radial basis functions, Eq.(7). This wave function produces the total energy E 
= -7.44859276608 a.u. for the ground 2 S— state of the °°Li atom. Only one electron spin-function 
Xi = a/3a — j3aa was used in these calculations. 



k 


m\{k) m,2(k) ra^(k) 








7fc 


1 








1 


0.146131429481911E+02 


0.416423958308045E+01 


0.390330393520923E+01 


0.754647967615054E+00 


2 


1 





1 


0.110681883475133E+03 


0.526120133598745E+01 


0.699723209902291E+01 


0.604874234309236E+00 


3 


1 


1 


1 


0.234004890152000E+03 


0.598950989816454E+01 


0.378246213278597E+01 


0.590931685390717E+00 


4 


2 





1 


0.116882861303726E+03 


0.413820964104499E+01 


0.322650638788059E+01 


0.590283203974117E+00 


5 











-0.112533836667930E+02 


0.404234943769060E+01 


0.402575246908481E+01 


0.136870242306376E+01 


6 








2 


0.124148253762594E+01 


0.772078877480103E+01 


0.208089611437233E+01 


0.906692873308709E+00 


7 








3 


0.129950444032051E+01 


0.125634586976915E+02 


0.167180833975207E+02 


0.111484879880926E+01 


8 


3 





1 


0.121940082563610E+02 


0.336154087416427E+01 


0.329089537444583E+01 


0.569461303170668E+00 


9 


2 


2 


1 


-0.815264312642982E+01 


0.315618697991587E+01 


0.293590042487235E+01 


0.642752252837086E+00 


10 








4 


0.231307852206996E-01 


0.216119418842614E+01 


0.354528921898964E+01 


0.848100637437686E+00 


11 


1 








-0.290323204459308E+02 


0.3031 11504960694E+01 


0.340407714984374E+01 


0.439831339492880E+00 


12 


1 





2 


-0.782805087195444E+01 


0.31987434232921 1E+01 


0.294880236447181E+01 


0.670201325893484E+00 


13 


4 





1 


0.466451849833029E+05 


0.257191421047538E+02 


0.198922117452300E+02 


0.745721864825035E+00 


14 


5 





1 


0.102312737866967E+03 


0.610735801244370E+01 


0.126738573321754E+02 


0.583192334409026E+00 


15 


1 


1 





0.496831723313492E+01 


0.214720891197913E+01 


0.139151014500059E+02 


0.332004403600092E+00 


16 


1 


1 


2 


0.331570408306929E+01 


0.313893628304852E+01 


0.288670468153774E+01 


0.656255394173358E+00 


17 


2 








-0.224488896165048E+03 


0.11 1494404746535E+02 


0.712253949382151E+01 


0.412242317869619E+00 


18 


2 





2 


-0.200426722868461E+01 


0.282336074847368E+01 


0.292360100848728E+01 


0.638644442755579E+00 


19 








5 


0.186252350189224E-02 


0.103900277658347E+02 


0.241603283627644E+01 


0.858660133717424E+00 


20 


2 


1 


1 


0.314011900712871E+01 


0.314575640023816E+01 


0.238785024725699E+01 


0.645792868876097E+00 


21 


3 





2 


0.164525876141841E+01 


0.318702139403640E+01 


0.348200049008609E+01 


0.663487519449685E+00 


22 


3 








0.549306690646211E+01 


0.313423005590041E+01 


0.376361954006700E+01 


0.494738425712773E+00 


23 


3 


1 


1 


0.489167893642778E+01 


0.310583281642480E+01 


0.299668864388672E+01 


0.617551331047668E+00 
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TABLE II: The probabilities (in %) of observing selected electron final states in the helium ion 
4 He + and tritium atom 3 H arising in the exothermic nuclear reaction Eq.(l) of the three-electron 
lithium-6 atom in its ground 2 S— state with slow neutrons. 



atom/state Is 


2s 


2p 


4 He+ 9.52386-lCr 2 


1.66605-10- 2 


6.60170-10- 4 


3 H 2.49810-KT 3 


3.05880-10- 4 


2.00251-10- 6 
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